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ABSTRACT 

We study anisotropic cosmology in Horava's gravity theory and 
obtain Kasner type solutions, valid for any number d of spa- 
tial dimensions. The corresponding exponents satisfy two rela- 
tions, one involving the marginal coupling A . Also, Horava's 
(super)renormalisable theory predicts (super)stiff matter whose 
equation of state is p = wp with w > 1 . We discuss briefly 
the implications of these results for the nature of cosmological 
collapse. 
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1. Introduction 



Recently, Horava has proposed [1] a candidate tlieory for gravity based 
on anisotropic scaling of space and time coordinates: 

x' ^ Ix' , t ^ Ft (1) 

where z is the scaling exponent. He has constructed an action for the metric 
fields invariant under the above scaling and also under foliation preserving 
diffeomorphic transformations. The action is required to have no more than 
two time derivatives. The kinetic part of the action is then universal and is 
characterised by a marginal coupling A . The potential part has numerous 
terms containing various powers and spatial derivatives of curvatures of the 
spatial metric. Hofava has invoked 'the principle of detailed balance' to con- 
strain such terms, but this seems unnecessary and may even be problematic. 
The action reduces to the Einstein action in the infrared (IR) if A ^ 1 
and z ^ 1 . The theory then has full space-time diffeomorphism symmetry 
and may, therefore, be a candidate for a renormalisable Einstein's theory of 
gravity. Hofava' s theory may also acquire an anisotropic Weyl symmetry at 
A = I where d is the number of spatial dimensions. 

Such a theory has many appealing properties. For example, it is ghost 
free since there are no more than two time derivatives. By construction, it is 
power counting renormalisable in the ultraviolet (UV) if z = d and is super 
renormahsable if z > d, so it is beheved to be UV complete. It contains 
many higher powers and derivatives of curvature, hence it may be able to 
resolve singularities. It singles out time, so the causal structure in UV is 
likely to be modified which may have non trivial implications for black hole 
physics. In such a theory, the speed of light generically diverges in the UV, 
so the horizon problem may perhaps be solved without requiring inflation. 
See [1, 2] for more details. 

Various aspects of such a theory are being studied actively. See, for 
example, [2] - [20]. In this letter, we focus on implications of such a theory in 
early universe cosmology where differences from Einstein's theory are likely to 
be manifest. Such imphcations have also been studied in [5, 6, 8, 9, 13, 14] for 
d — 3 homogeneous isotropic FRW universe. It is found that scale invariant 
cosmological perturbations can be generated without requiring inflation if 
the scale factor a{t) evolves as ~ with n > | ; it is also found that there 
can be a bounce in the early universe if the spatial curvature is non zero. 
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The scale invariance of perturbations is due to the modifications of dispersion 
relations arising from anisotropic scaling symmetry and, hence, is likely to be 
a generic feature of Hofava's theory. The bounce is due to non zero spatial 
curvature of the FRW universe and due to higher powers of curvature in the 
action. The bounce is thus a possible but non generic feature of Hofava's 
theory. For example, it is absent for spatially fiat universe. 

Our main motivation here is to find the implications of Hofava's theory 
which differ from those of Einstein's theory and which are not crucially de- 
pendent on the spatial curvature. We find that generic equation of state for 
the matter in the UV is ^ = w = ^ . This is independent of whether the 
spatial curvature is zero or non zero. Thus, in the early universe, w = 1 for 
renormalisable theory whereas w > 1 for super renormalisable theory. The 
corresponding matter is sometimes referred to as (super) stiff. 

We consider the general action given in [6] and study the evolution of 
a homogeneous anisotropic universe with d spatial dimensions. We obtain 
anisotropic Kasner type solution in the limit where universe is collapsing to 
zero size. The exponents in the corresponding scale factors satisfy two rela- 
tions, one of them involving A . The marginal coupling A may be different 
from 1 in the UV, and may even be close to | where the theory may acquire a 
Weyl symmetry. Such a behaviour of A and the presence of (super) stiff mat- 
ter have interesting implications for the nature of collapse which we explain 
briefly. 

This paper is organised as follows. In section 2 we present the set up. In 
section 3 we discuss the dispersion relation and the consequent equation of 

state. In section 4 we present the equations of motion and some solutions, 
and discuss briefly their implications. In section 5 we conclude with a brief 
summary and a few comments. 

2. Ansatz for action and metric 

In Hofava's theory, the fields are the lapse function A'", shift vector , 
and the spatial metric gij . The scaling dimensions of various quantities in 
momentum units are: 

[x']^-l , [t] = -z , [N] = [g,^]^0 , [N']^z-1. 

The action S = Sk + Sy is required to contain no more than two time 
derivatives, and to be invariant under the scaling in equation (1) and foli- 
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ation preserving diffeomorphism. The kinetic part Sk of the action is then 
universal and may be written as 

Sk = ^ J dt d^x [KijK'^ - XK^) (2) 

where is a parameter with dimension [k^] = z — d , X is a dimensionless 
parameter, the spatial indices i,j,--- = l,2,---,d are to be lowered or raised 
using Qij or its inverse g''^ , 

Kij = ^ {9^, - V,7V, - V.-iV.) , K = g'^ K,j , 

and the covariant derivatives, as well as curvature tensors below, are all with 
respect to gij . For z = d , the parameter k, becomes dimensionless and the 
theory is power counting renormalisable; for z > d it is super renormalisable 
[1]. Our interest is in the d — 3 case, but most of the expressions below are 
valid for any value of d . 

The potential part Sv of the action contains various powers and spatial 
derivatives of the Riemann tensor Rijki , equivalently of the Ricci tensor Rij 
in the d = 3 case. It suffices our purposes to write Sy symbolically as 

Sv= / dtd'xN^ L, + ei?+ECni?"+ E/^P^^v^^' (3) 

where the first sum denotes various powers of curvature tensor and the second 
denotes various derivatives acting on various powers of curvature tensor. The 
upper limits of {n, p, q) depend on the value of z . For the renormalisable 
case, for example, z = d and = z , + 2q^ + 2 = 2z . 

In [1], Hofava invokes 'the principle of detailed balance' which will con- 
strain the above form for Sy ■ For example, in o? = 3 case, Sy will not 
contain terms and the coefficients of various terms in Sy depend only on 
three new parameters. By construction, the resulting action is not the most 
general one. However, quantum corrections may not obey the principle of 
detailed balance and may induce other possible terms. This principle may 
even be problematic since the corresponding static spherically symmetric so- 
lutions reduce to the IR ones only on scales beyond the cosmological horizon, 
and not on smaller scales where Einstein's theory has been well tested. To 
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rectify this problem requires going beyond the detailed balance. See [5, 6, 7] 
and, in particular, [10] for detailed discussions of these issues. 

For these reasons, we will not invoke detailed balance in this paper, and 
consider the general form for the action Sy- The most general form of Sy for 
z = d = 3 is given, for example, in [6] where the corresponding equations of 
motion arc also obtained. These equations arc very long and, hence, are not 
presented here but will be used for the present cosmo logical study. Note that, 
for such a study, one can set N — 1 and N"^ — in the equations of motion 
with no loss of generality; and, also that these equations are applicable for 
any value of d as we later explain. 

Here, we consider only spatially curved, homogeneous, isotropic universe; 
or spatially flat, homogeneous, anisotropic universe. The line element of a 
spatially curved, homogeneous, isotropic universe may be written as 

ds^ = -dt"^ + d^lj^ (4) 

where a{t) is the scale factor, (iS^ ^ is the line element oi a d - dimensional 

space of constant curvature, and k — +1, —1, for positive, negative, or 
zero curvature. The Hubble parameter H is defined hy H — ^ where an 
overdot denotes time derivative. 

The line element of a spatially flat, homogeneous, anisotropic universe 
may be written as 

ds" = -de + Y,al {dx'f (5) 

where ai{t) are the scale factors. The corresponding Hubble parameters 
hi are defined by /ij = ^ . Also, define the geometric mean a of the scale 

factors by a'^ — Ylidi- Then H — ^ — ^ J^ihi is the average of the Hubble 

parameters hi . 

With the above definitions, the conservation equation for a matter source 
with pressure p and density p is given in both of the above cases by 

p + dH{p + p)^0. (6) 

If the equation of state is given hy p = w p then we have p = Po a"''^^'^"'^ 
where po is an initial value. 

3. (Super)stiff matter 
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Consider now matter sources, e.g. radiation, and their equations of state. 
The matter action which is invariant under the scahng in equation (1) will 
lead to a modified dispersion relation in the UV, typically of the form o;^ ~ 
k"^^ [1, 2]. From the standard statistical mechanical methods using such 
a dispersion relation, it follows that the dependence of free energy F on 
temperature T is of the form F ~ T^^~ [2, 17]. For renormalisable theories 
z = 1 in the IR and z = d in the UV. It then follows that F ~ T^+d- 
at low temperatures and F ~ 7"^+^ at high temperatures. As noted in 
[2], similar free energy behaviour at high temperature appears also in string 
theory. We further note here that similar free energy behaviour, at both low 
and high temperatures, appears also in the context of a particular version of 
generalised uncertainty principle [21]. 

With free energy F ~ T^+f , it follows upon using thermodynamical 
relations that the corresponding equation of state is given by p — w p where 
w — ^ . Thus for radiation in = 3 , we have z — 1 and w = ^ in IR. We 
have z = d in the UV for renormalisable theories which then implies that 
w = 1 , ^ the corresponding matter sometimes referred to as stiff matter. 
Also, z > d for super renormalisable theories which then implies that w can 
be > 1 , the corresponding matter sometimes referred to as superstiff matter. 

Such an UV dispersion relation, namely uP' ~ A;^^ , is ubiquitous in 
Horava's theory and arises from an underlying principle: it is a consequence 
of invariance under the anisotropic scaling in equation (1). Also, it is in- 
dependent of whether spatial curvature is zero or non zero. Thus, Horava's 
theory can be taken to predict that early universe, and more generally UV 
regime, is dominated by matter whose equation of state is given hy p — wp 
where w = ^ = 1 for renormalisable theories and > 1 for super renormalis- 
able theories. ^ 

4. Equations of motion, solutions and their implications 

Consider now equations of motion. They are given in [6] ior d — 3 . Con- 
sider, for any d , the general form of the contributions of various terms in 

^That to = 1 for radiation in the UV is also pointed out in [18] which appeared while 
writing up this paper. 

^The idea that early universe must be dominated by w > 1 matter also appears in 
different contexts. For example, see [24, 25] for the w = 1 case; see [26] and references 
therein for the w > 1 case. 
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the action to the equations of motion. For the cases of interest here, namely 
where the hne element is given by equation (4) or (5), we observe the follow- 
ing: 

(i) Matter source with an equation of state p = w p will contribute terms 
oc a~'^(^+"') in the equations of motion, see equation (6). 

(ii) Consider terms of the form i?" in Sy ■ For the isotropic case, it is easy 
to see that they contribute a term oc /c" a~^" in the equations of motion. It 
thus follows that such terms act as sources with equations of state p = w p 
where to = ^ — 1 and p = Cnk" ci~^^ . The constant (7„ depends on the 
index structure of i?" terms and their coefficients in Sy ■ For the spatially 
flat case, ^ = and the corresponding contributions all vanish. 

Note that n = 1 for i? term and w — ^ — ^ — — | for (i = 3; n = 2 for 
R"^ term and -u; = | — 1 = | for d = 3; and, formally, n = for cosmological 
constant term and w = —1 for any value of d . For the term i?"* with 
highest power of curvature, we have = z = d for the renormalisable case 
and w — ^ — 1 — 1 for any value of d . 

(iii) Rijki for a constant curvature space is given in terms of gij , and the 
scale factor which depends on t only. Hence covariant derivatives acting on 
curvature tensors will all vanish. Therefore the terms in the second sum in 
equation (3) do not contribute to the equations of motion. 

(iv) The kinetic part Sk of the action is universal for any values of d and 
z . Hence, the corresponding terms in the equations of motion are just those 
given in [6]. 

Using the observations (i) - (iv) above and the expressions given in [6], 
we can now write the equations of motion. 

Isotropic case 

For the isotropic case, the metric is given in equation (4) and the equations 
of motion may be written as 

d(Xd-l) = 2^ Y.P (7) 
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(Xd-l) {H + dH') = Y.^P-p) (8) 



where H — ^ and H = ^ — . The sum J2 in the equations above 
denotes contributions from matter source, and also those from i?" terms in 
Sv for which p = — 1^ p and p ~ k'^ a~^" where the constant C„ 
depends on the index structure of i?" terms and their coefficients in Sy ■ 
Sec [6] for exphcit expressions for the d = 3 case. 

We will assume in this paper that Xd > 1 , since this is the case in 
Einstein's theory for which A = 1 , and also that A is not arbitrarily close to 
I . Replacing by (^^jEr) ^ then renders equations (7) and (8) identical 
to those in Einstein's theory. 

The evolution of the scale factor is then straightforward to understand. 
Let a — > (cxd) in the limit t — > (cx)) . The evolution is then dictated by 
those sources for which w is largest (smallest). If the total coefficient of the 

2 

dominant sources is positive then it follows that a{\) ~ . 

If the spatial curvature is non zero then sources arising from curvature 
terms can have negative coefficients. Then the total coefficient of the domi- 
nant sources can be negative and this will generically lead to a bounce in the 
evolution of a{t) . The details, and even the presence itself, of the bounce 
depend on the nature and strength of other sources present and can only be 
obtained by further analysis incorporating these data. 

Let z = d = ?) . Horava's theory then predicts the existence of stiff 
matter in the UV for which w = \ . Assume the total coefficient of the 
sources with tt; = 1 , which are the dominant ones in the limit a — > , to 
be positive. This will be the case for spatially fiat universe for which there 
are no contributions from curvature terms. It then follows that there is no 
bounce and that a(t) ~ t?- in the limit t ^ . 

We now make a remark. It has been shown in [6, 8, 13, 14] that, in 
Horava's theory, scale invariant primordial perturbation spectrum can be 
generated in the UV with an additional scalar field and without requiring 
inflation. Scale invariance arises, essentially, from the dispersion relation for 
the scalar field in the UV which is of the form a;^ ~ A;® . For the desired 
dynamics of the perturbations thus generated, it is also required that H'^a^ 
be an increasing function of t , or equivalently that /°° ^ converge, which 
is taken to imply that the scale factor a evolve as ~ t^ with n > | . See 
[6, 8, 13, 14] for details. 
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However, as described above, it is likely that a{t) ~ in the UV. 
Although this violates the requirement n > ^ , there may be no adverse effect 
on scale invariance of the spectrum since H^a^ may still be an increasing 
function of t because of the presence of other sources in equation (7) which 
will become important as t increases. This is plausible but, nevertheless, it 
is desireablc to study in detail the effects of a{t) ~ ts in the UV on the 
scale invariance of the spectrum obtained in [6, 8, 13, 14] in Hofava's theory 
without requiring inflation. 

Anisotropic case 

For the spatially flat anisotropic case, the metric is given in equation (5). 
There are no contributions from i?" terms and the equations of motion may 
be written as 

X(fH^-J2ihi? = (9) 

i 

{Xd - 1) (hi + dHhi) = {p- p) (10) 

where hi = ^ , H = ^ J2ihi , and hi = ^ - (^Y . Note that summing 
equation (10) over i gives 

(Xd-l) (^H + dH^) ^ {p - p) . (11) 
We have H = ^ from the definition a'^ = Hi o^i • It then follows that ^ 

hi-H = Aia-'^ (12) 
d{Xd-l)H^ = 2K^p + A^a-'"^ (13) 

where Aj are initial values satisfying I]i A = and = Z]j(A)^ • Once the 
equation of state p{p) is given then, in principle, p{a) can be obtained from 
equation (6), a{t) from equation (13), and ai{t) from equation (12). Note 
that the initial values A encode anisotropic initial conditions, e.g. during a 
collapse, and also that Ai may be thought of as a source with equation of 
state p = w p where w = 1 and Po = ^ > see equations (6) and (13). 

^Equations (10) and (11) give equation (12). Using hi = dH gives the constraint 
= . Substituting hi in equation (9) then gives equation (13). 
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Consider now the dynamics of the evolution, assuming the equation of 
state to be p = wp . Equation (6) then implies that p = po a''^('^+'^) where 
Po > is an initial value. The evolution in the limit a — > oo , namely large 

2 

universe hmit, is similar to the standard one where a ~ t . The effect 
of A is unimportant unless A is arbitrarily close to ^ . 

Consider a universe collapsing to zero size, i.e. a— >0,ast^O. Let 
the scale factors ~ t"' in this limit. We study the following cases. 

w >1 

In the limit a ^ , 2k'^p ~ a"'^(^+"') > A^a"^'^ in equation (13) since 
w > 1 . It is then straightforward to show that 

2 ™— 1 2 

a ~ tda+^ , tti^Ci e" *^ ^^0+^ (14) 

w—l 

where q and c are constants. Thus, since t^+i — > in the limit t — > , it 
follows that the exponents a* in ~ t"' are all equal, and are independent 
of the initial values . Hence, the collapse is isotropic and stable under 
perturbations [26]. 

w <1 

In the limit a — > , the right hand side of equation (13) becomes B^a"'^^ 
where B^ ^ if w < 1 and B^ = 2kVo + if w = 1 . It is then 
straightforward to show that 

. . c.- = i-| (A-i)' . (15) 

This is a Kasner type solution. The exponents a* depend on initial values 
A'' and, since Z^j = , satisfy the relations 

Ea-^l, A-.i:("-)^4+(A-i)|. (16) 

For w < 1 we have B"^ = A^ , hence X = A is the only possible value. For 
w = 1 we have B'^ = 2«;^po + A^ , hence < ^ < 1 and ^ < X < A . 
Clearly X ^ A if 2k'^po < A^ . Also X ^ | if 2k'^po > which is the 
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only possibility in Einstein's theory, or if A — ^ which is a new possibility 
in Horava's theory and is valid for any values of po and . 

Kasner type solutions, in particular the value of X and the dependence of 
a* on initial values, provide an insight into the stability of the cosmological 
collapse process under generic curvature and/or anisotropic perturbations; 
namely, an insight into whether the collapse will be isotropic or anisotropic, 
whether smooth or will exhibit chaotic oscillatory behaviour, et cetera. 

Since cuj depend on the initial values Ai , the collapse will be generically 
anisotropic. Consider d = 3 case. The exponents cuj must satisfy the 
constraints in equation (16). U X = 1 then one of the ttj must be negative. 
Then, under curvature perturbations, the collapse will not be smooth and 
will exhibit chaotic behaviour. If X is sufficiently close to | then no 
can be negative and the collapse will be stable and non oscillatory under 
perturbations. 

In Einstein's theory A = 1 and, hence, smaller values of X may result 
only through smaller values of which necessarily requires stiff matter 

with w = 1 . Stability also results if superstiff matter with w > 1 is 
present as follows from equation (14). See [26] and the references therein for 
a thorough discussion of these issues. 

In Horava's theory, on the other hand, A is typically different from 1 in 
the UV. This theory may acquire an anisotropic Weyl symmetry if A = ^ , so 
it is possible that A — > ^ in the UV. If so then X may be naturally small even 
without stiff matter. But (super)stiff matter with w > 1 is also naturally 
present in this theory. This makes it more likely that collapse process is 
stable and non oscillatory. 

However, spatial curvature terms of high order are also allowed in Horava's 
theory. As explained in remark (ii) in section 4, this order is closely linked 
to the scaling exponent z, to which is also linked the presence of (super)stiff 
matter. Curvature terms typically lead to destabilising effects under curva- 
ture perturbations and our preliminary analysis indicates that they may be 
comparable to the stabilising effects of (super) stiff matter. However, there is 
also the stabilising effect that results if A is close to ^ in the UV but no com- 
parable, A— dependent, destabilising curvature effects seem to be present. It 
is therefore possible that the sum total of all these effects in Horava's theory 
results in a stable and non oscillatory collapse. Clearly, further analysis is 
necessary but it is complicated and quite involved, and is beyond the scope 
of the present work. 
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5. Conclusion 



We now summarise briefly and make a few comments. Our main motiva- 
tion for the present study is to find tlie implications of Hof ava's theory which 
differ from those of Einstein's theory and which are not crucially dependent 
on the spatial curvature. The implications of Hofava's theory we find which 
differ from those of Einstein's theory are: (1) The UV regime in the (su- 
per)renormaIisabIc case is dominated by (super)stiff matter, namely matter 
with w > 1 . (2) I]j(Q!*)^ , where a* are the Kasner exponents, can be dif- 
ferent from and smaller than 1 even without stiff matter present. These two 
implications are generic. (3) Equations of motion contain curvature terms of 
the form C„ ^" a~^" where n < dior the renormalisable case. The constants 
Cn can be positive or negative. The curvature terms may hence lead to a 
bounce in the evolution of the scale factor a{t) . The bounce is possible but 
non generic since it depends on the sign and the magnitude of C„ k"' as well 
as the nature and strength of other sources present. 

We now make a few comments. Hofava's theory is beheved to be UV 
complete. Also, it contains higher powers and derivatives of curvature. It 
is then reasonable to expect that this theory will also resolve the singular- 
ities. The big bang singularity is indeed absent if a{t) bounces back. But 
the bounce is not a generic feature. If there is no bounce then, as can be 
inferred from the present solutions, the big bang singularity is present. Note 
that black hole singularity is also present in all the static spherically sym- 
metric solutions to Horava's theory studied so far, see e.g. [7, 10, 11, 12, 15], 
although these solutions do differ from those in Einstein's theory. 

It is possible that, regarding the presence or absence singularities, this 
is the most one can see using classical action in Hofava's theory and that 
quantisation of the action is necessary to see any further. 

There is a similar situation in string/M theory. No classical string/M 
theory action so far has led to generic absence of singularities. However, it is 
likely that when the temperatures become comparable to string theory scales 
the classical description of the universe, and similarly of black holes, must 
be abandoned and string theoretic description should be used, see [27, 28]. 
But the details of such a description, or of how exactly the singularities get 
resolved, are not fully known at present. However, due to entropic reasons, 
the universe at this stage seems likely to be dominated by stiff matter for 
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which w = 1 [24, 25]. ^ 

In this context, note that Hofava's theory also predicts generically the 
presence of stiff matter in the early universe. The presence of such matter 
and/or the scaling arguments using equation (1) lead to the high temperature 
behaviour of the free energy F ~ [2]. This behaviour may be taken 
to signify that the spectral dimension of the spacetime in the UV is 1 + 
1 . This is shown to be the case for Hofava's theory in [2]. Such an UV 
spectral dimension is also observed in many candidate theories for quantum 
gravity e.g. causal dynamical triangulations, quantum Einstein gravity, spin 
foam theory, and string theory. This may perhaps be the case also for loop 
quantum gravity. It is indeed argued in [5] that two dimensional effective 
gravitational theories in the UV may be a generic feature of UV complete 
theories of gravity. See [5] for more discussions and an extensive list of 
references. 

If these similarities are more than just coincidences then it may be that 
there are also similarities in the ways in which singularities get resolved in 
any of these theories and in Hofava's theory. It is therefore important to 
study if singularities can be resolved by quantising the action in Hofava's 
theory. Such a resolution, besides being important on its own, may also 
provide insights into the theories mentioned above. 

Acknowledgement: We thank the referee for many helpful suggestions 
which improved the presentation of the paper. 
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